Structure of glassy lithium sulfate films sputtered in nitrogen (LISON): 
Insight from Raman spectroscopy and ab initio calculations 
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Raman spectra of thin solid electrolyte films obtained by sputtering a Li2S04 target in nitrogen 
plasma are measured and compared to ab initio electronic structure calculations for clusters com- 
posed of 28 atoms. Agreement between measured and calculated spectra is obtained when oxygen 
atoms are replaced by nitrogen atoms and when the nitrogen atoms form bonds with each other. 
This suggests that the incorporation of nitrogen during the sputtering process leads to structures 
in the film, which prevent crystallization of these thin film salt glasses. 

PACS numbers: 66.30.Dn,66.30.Hs 



I. INTRODUCTION 

Glassy thin film electrolytes are materials of consider- 
able technological interest. They are used in the design 
of modern solid state batteries, electrochemical sensors, 
supercapacitors, and electrochromic devices. Most of the 
materials nowadays are fabricated using the sol-gel meth- 
ods, but new possibilities are currently explored by us- 
ing sputtering techniques. An example are thin films 
produced by sputtering a 0.75Li2O-0.25P2O5 target in a 
nitrogen plasma (LIPON,^), which show ionic conductiv- 
ities of about 3 X 10~^ri~^cm~^. This material has been 
successfully used in microbatteries^ and electrochromic 
systems Through the mixed network former effect the 
ionic conductivity could be further increased to about 
9 X lO^'^^^^cm-i by taking a 0.75Li2O-0.25[0.2SiO2- 
O.8P2O5] target^ 

The sputtering technique is in particular interesting, 
since it extends the glass forming range. This has re- 
cently been observed in the ternary Li20-B203-Li2S04 
systems,^ where the sulfate content could be increased to 
amounts that would lead to crystallization when prepar- 
ing glass by melt quenching. In fact, the borate compo- 
nent can be fully eliminated by sputtering a Li2S04 tar- 
get in a nitrogen containing plasma, leading to an amor- 
phous material referred to as LISON. ^ With respect to 
the target material, one could speak of a "salt glass". 
It is also found that sputtering in argon plasma does 
not lead to an amorphous film, indicating that nitrogen 
plays a crucial role in the glass forming ability. Indeed, it 
was shown that nitrogen is partly incorporated into the 
material^ with a sulphur to nitrogen ratio of about 2:1. 
At present it is not clear how the structure of these salt 
glasses is built up, in particular in which way the incor- 
poration of nitrogen takes place and how this leads to a 
stabilization of the amorphous structure. 

The aim of this work is to get insight in the structural 
arrangements of this new glass system, how the nitrogen 
is incorporated in the structure and how this is related 
to the increased stability of the material towards crystal- 
lization. To this end we study LISON films by Raman 



spectroscopy and compare the experimental vibrational 
spectra with predictions from ab initio electronic struc- 
ture calculations based on density functional theory. Cal- 
culations have been performed for pure Li2S04 clusters 
and nitrogen doped Li2S04 clusters. 



II. EXPERIMENTAL 

The Raman experiments were performed at room tem- 
perature in argon atmosphere using a micro- Raman setup 
with an argon-krypton laser tuned to the 514 nm line as 
excitation source. With this experimental set-up we have 
a depth resolution of 14 /im. The power at the sample 
was 21 mW for the depolarized configuration (perpen- 
dicular direction of polarization of the incoming beam 
with respect to the polarizer in front of the detector) 
and 7 mW for polarized configuration. The integration 
time of a single measurement was 300 s. The analysis 
of the spectrum is based on the median of 27 spectra. 
Measurements on different spots of the sample gave the 
same spectrum, showing that the film was homogeneous. 

The LISON film had a thickness of 4 /j,m and was pre- 
pared by sputtering a Li2S04 target in nitrogen plasma 
at pressure 1 Pa;^ Due to the rather strong Raman re- 
sponse of the silicon substrate and the small thickness 
of the film, the measured spectra are dominated by the 
silicon signal. In order to obtain the response of the film, 
a reference spectrum from a clean silicon substrate was 
subtracted. In this procedure it was important to ensure 
that the substrate had the same orientation with respect 
to the incident polarization, wherefore we made use of the 
polarization dependence of the intensity of the 521 cm~^ 
Raman line of silicon. 



III. AB INITIO CALCULATIONS 

Electronic structure calculations were performed for 
Li2S04 clusters with 28 atoms (4 Li2S04 units). Start 
configurations were created based on the conception that 
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the lithium sulfate salt glass consists of intact SO4 units 
with lithium ions in between. Using the GaussView3.0 
program,— SO4 tetrahedra were constructed with random 
orientation. Lithium ions were added under the con- 
straint of minimum interatomic distance of 1.6 A. 

The starting configurations were geometry optimized 
using the Gaussian03 suite of quantum chemistry 
programs^*' (for the theoretical background underlying 
the methods see^^). First, the semi-empirical PM3 
method was applied, followed by a Hartree Fock calcula- 
tion with the 6-31G* basis set. These calculations took 
about 200 CPUh and ~ 200 iterations. For the final opti- 
mization using the hybrid density functional B3LYP and 
the same basis set only a few iterations were needed. 

In order to monitor the influence of nitrogen, modi- 
fied Li2S04 clusters were generated by replacing oxygen 
atoms by nitrogen atoms in the optimized Li2S04 cluster 
configurations. These modified clusters were again opti- 
mized by employing Hartree Fock and B3LYP with the 
6-31G* basis set. 

From the 2nd and 3rd order derivatives of the poten- 
tial energy surface (PES) at the minima found for the 
optimized configurations, the Raman spectra were calcu- 
lated. It is clear that the amorphous film will exhibit a 
large number of local structures corresponding to differ- 
ent minima of the PES. By optimizing different starting 
configurations it is possible to scan a part of the PES and 
find minima, whose associated clusters are good represen- 
tatives of the local structures of the LISON film. This 
is evaluated by comparing calculated vibrational spectra 
with the experimental. 



IV. EXPERIMENTAL RESULTS 

Figure [1] shows (a) the polarized and (h) depolarized 
Raman spectra for the LISON film. The bare silicon 
(dashed line) and the sample spectrum (film plus sub- 
strate, solid line) are shown in the insets of Figs. [T^,b. 
The bare silicon spectrum was scaled to fit the charac- 
teristics between 250 cm~^ and 350 cm^^. The main 
feature of the resulting difference spectra in Figs. [T^,b 
is a band between 900 cm~^ and 1035 cm^^. A second 
weaker band is found at 1225 cm~^ with a shoulder at 
slightly lower frequencies, around 1150 cm^^, in the po- 
larized spectrum in Fig. [T^. In addition two more bands 
may be discerned at 415 cm~^ and at 645 cm~^, in par- 
ticular in the case of the polarized spectrum, where the 
intense silicon mode at 521 cm~^ is absent and does not 
hamper the subtraction procedure. The intensity gap be- 
tween 440 cm~^ and 500 cm^^ in the depolarized spec- 
trum results from the fact that, in order to resolve the 
weak difference signals, the integration time had to be 
chosen sufficiently large. As a consequence the intensity 
of the silicon mode at 521 cm~^ present in the depolar- 
ized spectrum (see inset in Fig.lTJj) saturates the detector 
in this spectral region. 

In order to assign the bands originating from the LI- 
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FIG. 1: (a) Polarized and (b) depolarized difference Raman 
spectra of LISON films on a silicon substrate. The difference 
signal is obtained by subtracting the spectrum of a bare silicon 
substrate from the spectrum of the sample. Inset shows the 
raw data from the sample and the substrate (dashed lines 
silicon, solid lines LISON film sample). The spectra have 
been offset for clarity. 



SON film in the difference spectra in Fig. [T] we com- 
pare the results to the Raman modes found in crystalline 
Li2S04 marked as bars in the figure. The Raman 
spectrum of the crystal is dominated by the symmet- 
ric breathing mode of 804^ at 1014 cm^^, which either 
is very weak or found at lower frequencies in the glassy 
film. That this main band covers a range of frequen- 
cies where there are no corresponding modes from the 
crystal, suggests that the structure of the amorphous LI- 
SON film and the Li2S04 crystal are different also on 
the short range order scale. The modes at 1123 cm~^ 
and 1127 cm~^ in the crystal (appearing as a single bar 
in Figs. [T^,b due to their small separation), as well as 
the mode at 1194 cm~^ seem on the other hand to be 
shifted to higher frequencies and can be related to the 
bands around 1150 cm~^ and around 1225 cm~^ in the 
experimental spectrum. Similarly, the two bands in the 
film around 415 cm^^ and 645 cm"-'^ have likely their 
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origin in the pair of modes at 447 cm"^, 513 cm"-'^ and 
the pair of modes at 617 cm^^, 665 cm^^ in the crystal, 
respectively. 

V. THEORETICAL RESULTS AND 
COMPARISON WITH THE EXPERIMENT 

For the comparison of the experimental and calculated 
spectra we will exclusively use the polarized Raman spec- 
trum shown in Fig.[T^, as it is not disturbed by the silicon 
mode at 521 cm~^. To obtain a smooth spectrum from 
the calculated Raman lines, each line j with Raman shift 
ujj and intensity Ij is replaced by a Gaussian peak func- 
tion Ij exp[-(w - WjO^/A^] with A = 20 cm'^ 

In addition, a rescaling of the calculated Raman shifts 
is in general necessary to obtain good agreement with the 
experimental spectra, depending on the particular mate- 
rial and method used. Commonly used literature data 
for these rescaling factors exist for some gas and organic 
molecules. They have been determined by adjusting cal- 
culated to measured frequencies for well defined modes. 
In our case such literature data are not available and we 
have to determine the rescaling factor by other means. 
To do this we focus on the SO4 breathing mode of the 
crystal, since one should expect it to be present also in 
the clusters and the film (although, possibly, with re- 
duced intensity). The calculated SO4 breathing mode 
frequencies differ slightly from each other due to the dif- 
ferent local surroundings of the 804^ anions in the clus- 
ters. Taking the average of them and determining the 
quotient with the corresponding crystal mode frequency 
gives a stretching factor of 1.11. Using this factor yields 
theoretical spectra which appear to be shifted to higher 
frequencies. An optimal overlap of calculated and mea- 
sured spectra is obtained with a slightly reduced factor 
of 1.09. This minor correction should not be surprising 
as some shift of the breathing mode frequency can be 
expected when considering the crystal and the film. To 
summarize, for the comparison all calculated spectra are 
stretched by this factor of 1.09. 
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FIG. 2: Comparison of two calculated Raman spectra for 
Li2S04 clusters with the experimental Raman data. 



However, there are two important differences : (i) The 
experimental spectrum shows a broad shoulder towards 
lower frequencies from the main band at 1002 cm~^ down 
to 900 cm~^. As there exists not a single vibrational 
mode between 760 cm~^ and 960 cm^^, this shoulder is 
absent in the calculated spectra, (ii) In the calculated 
spectra there is intensity due to a mode at 1080 cm~^, 
while almost no intensity is seen around this frequency 
in the experimental spectrum. As a first hypothesis de- 
viation (ii) might be related to the finite size of the cal- 
culated cluster. A closer inspection of the eigenvector 
belonging to this mode reveals that it contains an oxy- 
gen atom, which moves in the normal direction of the 
surface of the cluster. Accordingly, this mode can be 
attributed to an artefact of the limited cluster size. 

Deviation (i) on the other hand might be due to struc- 
tural changes when nitrogen is incorporated in the struc- 
ture, which can be investigated by comparing experimen- 
tal spectra to calculated spectra from the nitrogen con- 
taining clusters . 



A. Li2S04 clusters 

Figure [2] shows the Raman spectra calculated from 
two optimized Li2S04 cluster configurations (dotted and 
dashed lines) in comparison with the measured Raman 
spectrum (solid line). 

As can be seen from the figure, the calculated spectra 
for the two different optimized clusters do not deviate 
much. This can be interpreted in the way that local 
arrangements participating in the dominating vibrations 
are still determined by the Li+ and 864^ ions and their 
mutual Coulomb interaction as in the crystal. As a con- 
sequence, the 804^ units are rather weakly coupled. 

In comparison to the measured Raman spectrum, 
the calculated ones resemble the overall band structure. 



B. Li2S04-N (LISON) clusters 

The substitution of two oxygen atoms by nitrogen 
atoms in a cluster has a strong effect on the calculated 
spectra, as can be seen from Fig. [31 where spectra of 
4 different clusters (labeled 1-4 in the figure) are com- 
pared to the measured Raman spectrum. This is not 
surprismg as haff of the 8O4" are replaced by 8O3N3-. 
The rather large amount of nitrogen correspond to re- 
sults from Rutherford backscattering experiments which 
yielded a ratio of about 2:1 of sulfur to nitrogen^^ and 
are also in agreement with previously reported valuesi^ 

A number of new bands appear in the calculated spec- 
tra from the nitrogen containing clusters. However, while 
the clusters without nitrogen lead to very similar spec- 
tra, the modified systems show pronounced differences. 
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FIG. 3: Comparison of four calculated Raman spectra (dotted 
lines) for LISON clusters with the experimental Raman data 
(solid lines). 



In general the spectra 2-4 do not compare well with the 
experimental data, while the calculated spectrum 1 is 
surprisingly close to the experimental data. In partic- 
ular, the broad shoulder of the main band is here well 
reproduced. The additional band at 1080 cm~^ is still 
present, but does not relate to nitrogen incorporation as 
discussed above. It is interesting to note that the opti- 
mal cluster 1 does not only compare most favorably with 
the experimental spectrum, but also has the lowest en- 
ergy (1.9 eV lower than the next lowest energy cluster 
associated with spectrum 3). 

Figure m shows the atomic configuration of the optimal 
cluster 1. In contrast to the clusters 2-4, it contains an in- 
teraction between two nitrogen atoms in close proximity 
at a distance of 1.47 A, which serves as a type of bridge 
between the two SO3N pseudo-tetrahedra. This suggests 
that the incorporation of nitrogen triggers the formation 
of structures, which, by hindering a reorientation of the 
anions relative to each other, prevent the crystallization 
of the film (or, more precisely, increases the free energy 
barrier for crystallization). 



VI. CONCLUSIONS 

By comparing experimental and calculated Raman 
spectra for LISON systems, we are able to show that 



a good representation of an amorphous LISON film is 
possible using small atomic clusters. Our investigations 
reveal that in the LISON structure nitrogen is replacing 
oxygen in the SO^" tetrahedra. In the configuration that 
results in Raman spectra in agreement with the experi- 
mental data we find a nitrogen-nitrogen bridge between 
two anions. This might be the underlying reason for the 
prevention of crystallization of the LISON material dur- 
ing the sputtering process, as these nitrogen bridges re- 
strain the reorientation of the anions and therefore sta- 
bilize an amorphous film structure. 




FIG. 4: Atomic configuration of the optimal LISON clus- 
ter (corresponding to spectrum 1 in Fig. |3}. Oxygen atoms 
are marked in red (dark gray, large spheres), sulfur atoms in 
yellow (light gray, mid-size spheres), nitrogen atoms in blue 
(dark grey, mid-size spheres), and Li atoms in silver (grey, 
small spheres). 
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